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The underlying mechanisms for one photon phase control are revealed through 
a master equation approach. Specifically, two mechanisms are identified, one 
operating on the laser time scale and the other on the time scale of the system-bath 
interaction. The effect of the secular and non- secular Markovian approximations 
are carefully examined. 



1 Introduction: 

The possibility of manipulating and controlling quantum systems by using quan- 
tum features is the goal of many modern control protocol ^ . An unusual type of 
control, one photon phase control (OPPC), has been the subject of recent atten- 
tion and discussion (see, e.g., Faraday Discussions 153^). In this case, control 
takes place by varying the relative phase of components of a weak pulse, while 
keeping the power spectrum of the light source fixed. A seminal proof ^ showed 
that one-photon phase control was not possible for isolated molecular systems in 
which control is over products in the continuum. A subsequent experiment on 
control of retinal isomerization in bacteriorhodopsin ^ motivated controversy 
and the need for clarification of conditions under which such control was possi- 
ble. This clarification, provided in Ref. 9 showed that control was possible for 
both closed systems and for open quantum systems under well defined conditions. 
Specifically, we showed the following for control of an observable O in a system 
defined by Hamiltonian 

(a) If the system is closed (i.e., not coupled to an environment), then one- 
photon phase control is possible if [^m, O] ^ 0. For example, control over iso- 
merization in an isolated molecule is possible since the probability of observing 
an isomer (e.g., cis or trans) is an observable that does not commute with ^m- 
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The situation is similar for control over product formation in radiationless transi- 
tions, such as intersystem crossing and internal conversion 

However, if [^m, O] = 0, then one-photon phase control is not possible in an 
isolated system. 

(b) If the system is coupled to an environment then, as above, control is pos- 
sible if [^M, O] 7^ 0. However, even in the case where [^m, O] = control may 
still be possible, in which case it is environmentally assisted. 

Reference 9 provided these general rules, derived in a somewhat different 
way below, but did not provide any insight into the way in which environmental 
assistance works, or the conditions under which such control is appreciable. This 
is the subject of this paper. Specifically, we demonstrate the significant role of 
non-Markovian dynamics and strong system-bath coupling in enhancing the ex- 
tent of one-photon phase control. To do so we utilize a master equation based 
approach to the open system case, which, as seen below, proves enlightening. 
The net result of this paper is the demonstration of the underlying mechanism 
responsible for environmentally assisted one-photon phase control. 

Note, also, as an aside, that we have also computationally demonstrated^^ 
one-photon phase control in a model retinal system, making a connection with 
the original experiments on related systems. 



2 One-photon phase control and closed quantum systems. 

Let us consider a molecular system M described by the Hamiltonian with 
Hm\^) =Ea\a). If the system state is described by the density operator pM, the 
expectation value of a given observable O of system M at time t can be calculated 
as 

{d{t))=Y,{a\d\b){b\pM\a). (1) 
If [^M, O] — 0, we simply get 

(^(0) = L(^I^I^)(^IpmI^)- (2) 

a 

For this case, only diagonal terms contribute to the time evolution of the expecta- 
tion value of the observable O. 

The essence of phase control is trying to manipulate the expectation of O by 
changing the phase of an applied laser pulse. When the laser pulse is taking into 
account, the Hamiltonian can be re- written as = +^^^ = — d • E(z, 0' 
with E(z,0 = ££(^,0 where 



'E[z,t) = j dcoe(co)exp -\(i){t--^ 



(3) 



is the amplitude of the field with polarization vector £, which can be expressed 
in terms of frequency spectrum e(co) = |e(co)| exp(i4)[co)] and, d is the transition- 
dipolar operator. We are interested in manipulating the time evolution of the 
observable O by means of varying the spectral phase (|)(co) only. 

In order to clearly understand the role of the environment in OPPC, it is il- 
lustrative to analyze in detail the circumstances under which ^ phase control is 
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not possible. The evolution of the density matrix, for unitary time-evolution, is 
governed by 

?ab = -iCO,,p,, + ^ X [HT{t)9cb - 9acHT{t)\ (4) 
^ c 

where p^^ = (a | pM | , hi^^ab = Ea-Eb and H^^ (t) = dab'Eiz.t) = {a\£-d\b)T.{z,t). 
Due the unitary character of the dynamics, we can integrate Eq. (4) out by means 
of a the following ansatz for the wavefunction 

withpM(o = i^(0)moi. 

Let us consider, e.g., that the molecule is initially = — oo) in the pure state 
l^"!), so we have b\ {t = — oo) = 1 and bk{t = — oo) = for A: 7^ 1. After formally 
integrating Eq. (5) for bm, we get 

bm{t) = f dsbi{s)cxp{mmAs)H^^{s). (6) 

If the perturbation is weak, for unitary evolution, one can assume that bi {t) = 1 
at all times, thus 



bm{t) = \dm\ I dcoe(co) I d^exp(i(co^,i -co)^), 
n J — 00 J — 00 



(7) 



where we have made use of Eq. (3) and defined e(co) = e(co)exp(icoz/<:). Since 
our interest is in the long time regime, we extend the upper integration limit to 
infinity and get 27r5(co^i — co) for the integral over s. Finally, we have 



/?^(oo) =27l^d^i|e(co^i)|exp 
n 



i ( (t)(co^i) + -^co^iz 



(8) 



which yields 



471 

Pab{t) = -2-d«id^i|e(co«i)||e(cO/,i)|exp(-ico^/,^) 
n 

X exp |i ^4>(co^i) - (t)(coz,i) + ^(co^i - co/,i) | . 



(9) 



From here, it is clear that no spectral phase information is present in the diag- 
onal terms, a = /?, and therefore it is impossible to manipulate the time evolution 
of the observable O in Eq. (2). 

Based on the equation (9), it is easy to track the vanishing of the spectral 
phase contribution to the fact that pM was constructed as the projector onto the 
state 1^(0)' i-^-' to the fact that the evolution is unitary. This can be understood as 
follows. Under unitary time evolution, the dynamics is generated by the unitary 
time-evolution operator {t) = 'rexp(— i fl^ dsH{s)/h), by means of the relation 

Pm(0 = ^ (OPm(0)L/^^(^). If p(0) represents the density operator of a pure state 
1^), then U{t) evolves |^) (in the direct Hilbert space) while U^t) evolves (^| 



This journal is ©The Royal Society of Chemistry [year] 



Foraday Discuss., [year], [yfo\], 1-11 |3 



(in the dual Hilbert space). Since the generator in the dual space, U"^ (t), is just the 
time-reversed version of the generator in the direct space, U{t), then the evolution 
in the two spaces interfere in the density operator in such a way that no phase 
information is encoded in the diagonal elements [see Eq. (9)] . This will no longer 
be true in the case of non-unitary dynamics (e.g., for open quantum systems) 
because of the time- symmetry breaking. This fact can be identified as a first 
contribution from the environment to OPPC (see below). 

In what follows, we explore in detail the possible ways in which the environ- 
ment can support the generation of OPPC. 

3 One-photon phase control and open quantum systems. 

The study of open-quantum-systems dynamics is a demanding and intricate task 
due to the various time/energy and coupling strength scales involved. In the case 
of molecular systems, it is common that one component of the system plays the 
role of the environment. For example, if one is interested in the electronic dy- 
namics, then the vibrational and rotational modes of the same system comprise 
the environment. In this situation, one expects that the coupling to the environ- 
ment is strong and the dynamics are non-local in time, i.e., it depends on the 
previous conformational states of the molecular system (non-Markovian dynam- 
ics) This non-local character of the dynamics has been proved to be vital in 
physical systems under the influence of external fields ^^'^^ Additionally, there is 
no apparent justification for neglecting the initial correlations between the part 
we name the system and the part we name the environment. All these features 
prevent us from applying standard approximation schemes and somehow hide the 
underlying physical processes. 

In order to identify the contribution from the environment to the one-photon 
phase control, we consider the following situation: a molecular system M is at 
equilibrium with its environment, a stationary/incoherent situation that is sub- 
sequently irradiated by a coherent laser pulse. In the first part of our analysis. 
Sec. (3.1), we neglect the initial correlation between the system and the envi- 
ronment and discuss how the non-unitary dynamics induced by the environment 
could assist OPPC. In the second part. Sec. (3.2), we include the initial corre- 
lations with the environment and show how they could enhance the role of the 
environment in assisting OPPC. 

3.1 Influence of the non-unitary dynamics 

In order to consider the influence of the environment on the system in a consistent 
way, we start by considering the total dynamics of the two parts by means of the 
total Hamiltonian 

H = Hm^H^^^He^ H^^, (10) 

where denotes the interaction between the system and the environment and 
characterizes the environment. Let us consider that the interaction between 
the system and the bath is of the general form 

H'^^ = Y,Ku^^u. (11) 

u 
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where Ku and are operators of the system and the environment, respectively. 
Note that Ku and could be non-Hermitian as long as is Hermitian. At 
this point we assume that the initial density operator pME of the system+bath 
factorizes at r = into a product of the system part pM and the environment part 
pE, that is as Pme(^o) = Pm(^o) ® Pe(^o). 

For this case, the non-Markovian evolution of the system density matrix ele- 
ments can be cast in the form (cf. Chapter 3 in Ref. 15) 

Pab = -'^^ab9ab ^ T^Y^ {^^^9cb - PacH^^) 

c 

t-to 

-^1 dx{M,^,^,(-x)ei^^^V,,(^-x) + M,,,,^(x)ei^^^V^,(^-x) d^) 



^dbA-^y""''" + M^,,,,(x)e^^^«^] pcd{t - X)} 



where M^^ ^^(0 determines the time span for correlations and are called the mem- 
ory matrix elements (cf. Chapter 3 in Ref. 15). They are defined by 

MabMt)=LCuAt)Ki"t^Kl'J, (13) 

and satisfy the relation M*^^j(r) = Mjc,^a(— 0- The reservoir correlation func- 
tion Cuv{t) is given by 

Cuvit) = i(<^>.(0^v(0))R-;^(<l>.)R(<l>v)R, (14) 

n n 

where {<t>u{t)^v{0))R = tr[pE,equ^M(0^v(0)]. The first line in Eq. (12) describes 
the unitary evolution and is completely equivalent to Eq. (4). Hence, if no phase 
control can be achieved from Eq. (4), no phase control can be achieved from this 
part of the evolution in Eq. (12). The remaining lines describe the non-unitary 
non-Markovian dynamics, characterized by the non-local evolution implied by 
the integration from to ^ — ^o- The most relevant feature of this second part 
of Eq. (12) is the fact that the off-diagonal elements affect the dynamics of the 
diagonal terms, despite the weak field condition. 

Consider then the case where the system is initially in equilibrium with the 
bath and the excitation follows. If we write Eq. (12) as pM = — iXpM — i^pu — 
iAaserpM, then our initial state is an eigenstate of the super operator 'rexp[—i fl^ ds{L{s) - 

ii>{s))], i.e., 'rexp[— i fl^ds{L{s) — i(Z>(^))]pM = A-pj^pM- Here L comprises the 
contribution from the system Hamiltonian Hm while i> provides for the effect of 
the bath on the system and Xiaser accounts for the action of the laser pulse. For 
high enough temperature and weak coupling to the bath [cf. Ref. 16 and refer- 
ence therein], one can show that in the long time regime, and in the absence of 
the laser, the density operator of the system is (where \n) are eigenstates of Hu)- 

PM,can = limpM(0 = ^6-^^(3 = ^ ^e-^'^P |^) (^| , (15) 

Here Z is the partition function Z = tr(e~^MP) ^nd |3 = 1/^b^, is an eigenstate of 
non-driven dynamics with eigenvalue 1, i.e., exp[— i fl^ d5'(X(^) — i(Z)(^))]pM,can = 
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-i{L- 


i2))pM,can, 






-i(2- 




-iAaserpM (0> 




-i(Z- 


i2))pil^(0 


-i£laserpM^(0- 



pM,can or equivalently, pM,can is a zero mode of the generator L{s) — i(b{s). Thus, 
in this description the laser excitation finds the system in an incoherent superpo- 
sition of energy eigenstates. 

For this initial state, under the weak field assumption, we can solve Eq. (12) 
perturbatively in the laser term, 

(16) 
(17) 
(18) 

The zeroth-order term p^{t) describes the dynamics in the absence of the laser 
and based on the description above, can be trivially obtained as (t) = pM,can. 
The first-order term {t) accounts for the dynamics induced by the laser pulse 
on the incoherent superposition of eigenstates pM,can. Since — i(X — ii>)p^^'> {t) = 
— i(X — i4^)pM,can = 0, the system-bath contribution is absent here. 

In deriving Eq. (9), we had assumed that Pm(0) = 1 1) (1 1, whereas in our case 
we have an incoherent superposition of eigenstates (15). Hence, by virtue of the 
linearity of the master equation (12), we can infer directly that 

Pabi^) = ^^l^e-^'^Pd«^d^Je(co«^)||e(co^„)|exp(-ico^^O 

^ ^ ^ (19) 

X exp |i ^^{(Jdan) - (^(CO^^) + ^(cO^n - (^bn)] } • 

Therefore, it is clear that from Eq. (19) that no phase information is encoded in 
the populations, a = b,of p^{t). 

The second-order term pj^^ (t) accounts for the evolution of the state prepared 
by the laser excitation in the presence of the system-bath dynamics. For the 
populations of p-^^ we have 

a(2) _ i V /^^/MR^(l) p,(1)z:.MrA 
Paa - ^ 2^ [^ac Pca - Pac ^ca J 

t-tQ 

Explicitly integrating Eq. (20) requires that we specify the system-bath interac- 
tions [Eq. (11)] and the bath nature. For real environments this integration is 
not possible and one has to appeal to specific models; in the Appendix we dis- 
cuss the most relevant features of the commonly used bosonic bath. Despite the 
sheer complexity of integrating Eq. (20), one can still extract some of the rele- 
vant physical contributions directly from the master equation (as an example, see 
Ref. 17 for a specific implementation). By inserting Eq. (19) into the right-hand 
side of Eq. (20), it is clear that phase control of the population will result from 
the system-bath interaction, which couples p^^^ to p^aa- This phase control is of 
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second order in both the radiation field and in the system-bath coupHng strength. 
Note that the radiation field (the first two terms) also couples p^^^ to p^aJ - This 
contribution, however, is third order in the radiation field and thus is negligible 
in our current considerations of second order effects. 

The time dependence of the M matrix elements signifies the non-Markovian 
character of the dynamics, which can be associated with the environment self- 
correlations. These self-correlations allows for a back and forth information flow 
between the system and the bath^^. If the dynamics is Markovian (see next case 
below), the information flow is unidirectional: from the system to the bath. This 
information flow can be understood as an entropy flux and in the case of Marko- 
vian dynamics results in an always increasing entropy of the system M. In the 
presence of external driving fields, non-Markovian terms are essential In our 
case, this back and forth flow allows for an efficient transfer of phase informa- 
tion via the bath. Additionally, the non-Markovian dynamics is characterized by 
slower correlations decay-rate allowing thus the presence of phase control for 
longer times 

3.1.1 Markovian Approximation — If the environment quickly "forgets" any 
internal self-correlations during its interaction with the system, that is, if any dy- 
namically established quantum correlations between parts of the environment are 
destroyed on a time scale much shorter than the characteristic time scale of the 
interaction-picture reduced density operator of the system, the Markovian approx- 
imation can be invoked. Mathematically, this means that the correlations Qv(^). 
and hence the M matrix elements, die off quickly as x > 0, and thus the x inte- 
gral in Eq. (20) picks up contribution only from p^^J{t), rendering the dynamics 
time-local. Moreover, it does not make a difference if we extend the upper limit 
of integration to infinity, in which case the matrix elements Tab^cd^ defined below, 
become constant. As a result, a Markovian master equation is obtained: 

Paa — \^ac 9ca " Pac J^ca J 

c,d ^ 

where 

rabM^) = ^ rdxe'^MabM^) = 3iY,Ki"^Kld rdxe'^CuA^). (22) 

JO JO 

This master equation retains the essential feature of the coupling of p^^J to pi^^ as 
in Eq. (20). As a result, even given the Markovian approximation, the population 
will show phase dependence of second order in the radiation field and of second 
order in the system-bath coupling. The Markovian approximation, which results 
in the time locality of the master equation, does not undermine the essential fea- 
ture of phase dependence in the population. However, Markovian dynamics re- 
sults, in general, in a stronger correlations decay-rate and therefore, the phase 
information will rapidly flow back to the bath. Hence, phase control may not he 
noticeable when the bath relaxation is very fast. 
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3.1.2 Markovian Secular Approximation — If, in addition, one now applies 
the secular approximation in Eq. (21), we get 

a(2) _ i V frrURjl) A'^)TrMR\ 

Paa — J^L^ y^ac 9ca " 9ac J^ca J 

(23) 

- 2^ [racM^ac)paJ {t) - rcaM^ca)pcc {t) . 

c 

This approximation decouples the evolution of the diagonal elements from that of 
the off-diagonal terms via the bath. Hence no phase control in the population can 
be observed, due to the artificial truncation of the off-diagonal terms. Thus, the 
secular approximation is not appropriate for the study of environmentally assisted 
phase control. Alternatively, from a physics perspective, if the secular approxi- 
mation is indeed valid for some system of interest, then our analysis shows that 
OPPC is not achievable in such a system. 

3.2 The significant infiuence of the initial system-bath correlations 

The presence of initial correlations between the system and bath, Pme(^o) 7^ 
Pm(^o) ^ Pe(^o). introduces an important new contribution in the master equation 
[Eq. (12)]. For completeness we write the entire expression. 



Pab = -iO^abPab + ^ {H^^Pcb - PacH^ 



-■^^(a|trE([^^,<i>^„pME(^0)] - [Ku^u, pE {to) pM {to)]) \b) 

(24) 

where Pe(^o) = trspME(^o) and Pm(^o) = trEpME(^o)- The new term in the dynam- 
ics exposes the initial correlations between the system and the environment. Note 
that this is an extension beyond standard master equations, which do not usually 
include the state of the environment explicitly. 

If Pme(^o) factorizes into system and bath terms, then it is clear that the terms 
in the second line in Eq. (24), the ones which account for the initial correla- 
tions, vanish. In the opposite case, this would imply that the initial stationary 
state should be re-derived taking into account this new driving term, i.e., the new 
stationary state will not be diagonal in the energy eigenbasis. The immediate 
consequence of this fact is that the initial state, before the coherent excitation 
takes place, will contain stationary coherences i.e. non-zero time independent 
off-diagonal elements of the system density matrix in the energy representation, 
which will allow for phase control. 

As in the previous case, we can write Eq. (24) as pM = —iLpu — iA)pM — 
^pM — iAaserpM and assume that the initial state is an eigenstate of the superop- 
erator 'T exp[-i fl^ ds{L{s) + Lo{s) - i'D{s))]. Here X, (D, Xi aser are defined as 
before and Lo denotes the contribution from the initial correlations given in the 
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second line of Eq. (24). Due to the initial correlations between the system and the 
bath, it is clear that in the eigenstates of the super operator 'rexp[—i fl^ ds{L{s) + 

Lo{s) — ii){s))] are no longer diagonal in the eigenbasis of ^m- For this case, 




where Z' = tr(e~(^M+^^^+^E)Pj ^j^^ partition function. Thus, in this case the 
laser excitation finds the system with stationary coherences of energy eigenstates, 
a concept that we have discussed in Ref. 14. The presence of these off-diagonal 
elements guarantees that the phases of the laser pulse contribute to the dynamics 
of the populations. 

Following a perturbative approach similar to that in Sec. 3.1, we can study 
the dynamics by means of the sequence 

it) = -i{L + io - i^^)pM,eq, (26) 

(0 = -i(x + io - im^M^ (0 - i wii^ (0, (28) 

Since the off-diagonal elements in the system energy eigenbasis are present from 
the beginning, phase control of the population will take place at the time scale 
of the radiation field. In other words, the phase information can "flow" directly 
into the population via the radiation, without having to wait for the system-bath 
interaction to first establish phase dependence in the off-diagonal elements. 

4 Concluding Remarks 

Environmentally assisted one-photon phase control arises in open systems when 
control is achieved over system properties that commute with the system Hamil- 
tonian. In such cases, phase control is not possible unless the environment partic- 
ipates in the dynamics. In this paper we have analyzed the origins of such control 
and shown that the effects of molecule-bath coupling on phase control are two- 
fold and take place on two different time scales. First, the thermalization of the 
coupled system-bath gives rise to non-zero off-diagonal elements in the initial 
molecular state, which can result in a direct phase control of the population; this 
control takes place at the timescale of the radiation field. Second, through the 
molecule-bath coupled dynamics, as modeled by Eq. (20) or Eq. (21), the phase 
information contained in the off-diagonals can "flow" through the bath into the di- 
agonal elements, resulting in an indirect phase control; this control takes place on 
the timescale of the molecule-bath interaction. The different timescales provide 
useful indicators to distinguish the two kinds of phase control in experiments. 

One final note is in order. Environmentally assisted one-photon phase control 
depends on the system-bath interaction. In the case where the bath is in equilib- 
rium this interaction will cause, long after the pulse is over, the system to relax 
to a laser-phase-independent equilibrium state as well. As such, and under these 
circumstances, environmentally assisted one-photon phase control will only sur- 
vive for a well defined time after the laser excitation is over. Rates and time 
scales over which this control persists is dependent upon, and requires numerical 
studies of, individual systems of interest. 
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Appendix: Harmonic Oscillator Bath 

Since the explicit calculation of the correlation functions C^v in Eq. (14) is not feasible for real baths 
such as solvents, one has to appeal to specific models of the bath and the system-bath interaction. 
The most common approximation is the normal mode approximation. In this case one assumes that, 
e.g., the vibrational modes can be described by small oscillations around the equilibrium point, thus 
allowing the use of the harmonic approximation. In this case, we replace the interaction term in 
Eq.dDby 

h"^^ = K{s)^Z) ^ K{s)Y,hy^Z^, (29) 

where s denotes the system coordinates, denote the system-environment coupling constants and 
Z = {Z^} denotes the environment normal mode coordinates. For simplicity, we have assumed that 
contains a single term and have dropped the index u. Note that since (Z^)r = 0, the expression 
for Cuv in Eq. (14) reduces to the two-point correlation function of the bath operator. 
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Thus, we have for Qv(0 = C(t) 

h(0 



C{t) = dco ^cos(coOcoth y2j~fj -isin(coO j coV(co), (30) 
where /(co) is the spectral density, 

Some typical spectral densities are, cf. Ref. 12, co^/(co) = 0(co) joco^e"^/^'^ or co^/(co) = 6(co) J^^^ ■ 
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